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Abstract 

We study the (zero temperature) quantum piezoelectric response of Harper-like 
models with broken inversion symmetry. The charge transport in these models 
is related to topological invariants (Chern numbers). We show that there are ar- 
bitrarily small periodic modulations of the atomic positions that lead to nonzero 
charge transport for the electrons. 

The Harper model can be interpreted as a tight-binding quantum Hamilto- 
nian describing the dynamics of non-interacting electrons on a two dimensional 
lattice in the presence of magnetic fields. It is known to have interesting Hall 
transport properties. Here we study the electric response of Harper-like mod- 
els to adiabatic changes in the hopping amplitudes. Changes in the hopping 
amplitudes have a natural interpretation as elastic deformation of the underly- 
ing lattice. As we shall show, such deformations can drive electron transport. 
We shall refer to this kind of response as piezoelectricity. Like the Hall con- 
ductance in the integer Hall effect |j, and in quasi-one dimensional systems 
[|, the Thouless pump |, §, the Ma gnus force Q, adiabatic charge transport 
in networks 0, adiabatic spin transport and adiabatic viscosity it is a 
transport phenomenon related to the adiabatic curvature and Berry's phases 

Let us first summarize the central findings: 1. Harper- like models with bro- 
ken time reversal and broken inversion symmetry have, in general, nontrivial 
piezoelectric response. 2. Appropriate periodic modulations of the atomic po- 
sitions lead to integral charge transport given by appropriate Chern integers. 
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This implies that an ac driving has a response with a dc component. 3. There 
are arbitrarily small periodic deformations that transport integral (and nonzero) 
charges over macroscopic distances. These periodic cycles trap level crossings 
in parameter space. 

The topological significance of piezoelectricity was noted in |^ , but vanished 
for the models considered there (Landau Hamiltonians). Thouless Q and sub- 
sequently Niu Q constructed one dimensional models of charge pumps where 
integral charge transport (given by a Chern number) is driven by one periodic 
potential sliding past another. Here, the model, the driving and the details are 
different, but the spirit is the same. 

We shall focus on a family of Harper models, H{t, fc, 77, (jj), which arises from 
tight binding models associated with a two dimensional triangular lattice. Each 
site of the lattice has a coordination number six and the basic plaquettes are 
triangles. Each up-triangle in the lattice is surrounded by three down-neighbors 
and vice versa. The magnetic flux through the up-triangles is 0/2-1-77 ^-^d 0/2 — 77 
through the down-triangles, (j) — 271 p/q, with p, q relative primes, and I77I < 7r/2 
is a measure of the asymmetry in the fluxes through the up/down triangles in 
units where the quantum of flux is 27r. The hopping amplitudes associated with 
the three basic vectors of the triangular lattice are tj ^ H, j — 1, 2, 3. The 
corresponding Harper model is: 

(^H{t, fc, 77, 0)*) (77) - (ii + hyne"f'/^) x^{n + 1) + 

2t2 COS(770 + fc2)*(77) + (il + tsynC'^/^^ X*(77 - 1). (l) 

X = exp(iA;i), ?/„ — ex.pi{n(l) + k2 — rj), ^'(71 + q) = '^{n) G C and x,yn are 
the complex conjugates of x,yn. k are Bloch momenta with ranges |fci| < Tr/q, 
1^2 1 < TT. The model was introduced in who studied the Hofstadter spectrum 
in the case t = (1, 1, 1). 

The class of models in Eq. (|l|) is the simplest among Harper-like models with 
interesting piezoelectric response. The simpler versions of the Harper model 
and, in particular, the classical Harper model on the rectangular lattice and 
its generalizations do not have interesting adiabatic piezoelectric response. 
The reason for this is that inversion symmetry needs to be broken. This is a 
fact about piezoelectricity that goes back to the brothers Curie Inversion 
symmetry is broken if 77 7^ mod tt. Inversion symmetry is preserved in the 
classical Harper model and the generalizations studied in ||l2| . 

Let \ip{t, k, rj)) be a normalized Bloch state of the Harper model in Eq. (|l]). 
Consider a closed loop 7 C IR^ in the space of hopping amplitudes. When 7 is 
traversed adiabatically (this, of course, subsumes that the gap remains open), 
the charge Q(7, fe, 77) transported from —00 to 00 in the ki direction for fixed 
k2 channel for each full band is given by Q 

Q(7,fe,77) = -Im / dk, {-^\Vt^p)-dt. (2) 

71" J-TT/q Jj Okl 
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The charge, if weU defined, is an integer — a Chern number. The total charge 
transported by the system is the sum over the relevant k2 channels and the 
occupied bands. When the system is an infinite two dimensional crystal then 
all the ^2 channels are relevant. On the other hand, for a strip of finite width 
with (possibly twisted) periodic boundary conditions, only a discrete set of 
values of k2 contributes. For reasons that shall become clear later, finite strips 
are the more interesting case. 

The difficulties in studying Chern numbers of model Hamiltonians ^, ^ Q 
(and this one is no exception) are: First, one needs to establish that the Chern 
numbers are well defined. For the problem at hand, this means that one needs 
to isolate a range of parameters fc2 , rj and for which the gaps surrounding an 
energy band remain open when t and ki run over their full range. Second, the 
Chern number may be well defined but zero, a case that is not very interesting 
for transport. For this not to be the case, the surface of integration in Eq. (||) 
must be protected against contraction. For certain transport properties such as 
those considered, e.g., in 0, 0, ||, the surface of integration had such a protection 
built in. This is not the case here. The cycle of deformations, 7, is a closed 
orbit in the three dimensional space of deformations, and such an orbit can be 
contracted to a point. If during this contraction the integrand in Eq. ^ remains 
continuous, the Chern number is zero. So, for the Chern number to be nonzero, 
the orbit of deformations 7 must trap level crossing. Finally, one needs to worry 
about global questions: Q(7, ^2,77) must be well defined for all of the relevant 
/c2 channels and must not sum up to zero. It turns out that the Harper model 
is rich enough so that everything actually happens there; there are good orbits 
and parameters where one finds nonzero quantized transport, but also bad ones 
where various bad things happen. 

The Bloch Hamiltonian, Eq. (^, is a homogeneous function of t of order 
one: iJ(r, k, Tj, 4)) = 1^1 H{i, k, Tj, (j)). The eigenvectors are independent of and 
the length of t therefore does not contribute to Eq. (||). We shall henceforth 
take t to be on the unit sphere. The three additional continuous parameters, 
k and rj are angular variables. Eq. (^ depends on five continuous parameters, 
{t, k, rj). The five dimensional parameter space is topologically the product of a 
two-sphere and a three-torus. 

To get one's hands on the Chern numbers for this model, one needs, as 
we have seen, to have good control over level crossings. One can use symme- 
try considerations to reduce the study of crossings from the full range of the 
parameters to a part of the parameter space. Indeed, there are three linear 
transformations of the parameters which are implemented by either unitary or 
anti-unitary transformations. These are 

{kj ^ kj + 2n/q}, {rj —f —rj, k^-k}, 
{fci ^ r, - fci - fc2 + (1 - {~m 0/4, ii ^ ts} (3) 

As a consequence of this, the spectral analysis of H{t, k, rj, 0), can be restricted 
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to the range: 



-'K/q<kj Kir/q, 0<?7<7r/2, h<tx. (4) 

We shall take ti and t2 to be our coordinates on the sphere of deformations. 
For the sake of concreteness we restrict ourselves to the positive quadrant tj > 
0, j = 1, 2, 3, and to the ground state of Eq. (|l|). We shall call the point on the 
unit sphere with tj = 1 "the j-th pole." 

Let r be the set of points where the lowest eigenvalue of Eq. (|l|) is degenerate. 
By r(fc2 = c) we shall denote the restriction of F to the subspace with fixed 
channel k2 = c and by F(fc2 = 0,7] — d) we denote the restriction to a fixed 
channel and asymmetry, etc. Recall that the von Neumann- Wigner rule 
says that a complex Hermitian matrix which depends on d parameters has, 
generically, eigenvalue crossings on a surface of dimensions d — 3. One therefore 
expects F to be two dimensional surfaces, F(fc2 — c) to be one dimensional 
curves and F(fc2 = c, 77 = d) to be isolated points. We shall see that this is a 
good guide to the behavior of the set of level crossings away from special points, 
e.g., the poles. For a generic point oit, the von Neumann- Wigner rule says that 
V{t) is a discrete set of points in k<E>ri space. At the poles we shall find, instead, 
that T{tj = 1) is a two dimensional surface. Of course, the poles are special 
points, and the failure of von Neumann- Wigner there is no source of concern. 

At the poles Eq. (|l|) can be diagonalized by hand. At the 2-pole the Hamil- 
tonian is already in a diagonal form. At the 1-pole it is diagonalized by plane 
waves and at the 3-pole, by plane waves up to an appropriate gauge transforma- 
tion. The restrictions of F to the poles, T(tj = 1), can be determined explicitly. 
More precisely, T{tj = 1) is the 2D set of points that obey: 

fci - (l + (-l)'')0/4, forj-1; 

fc2 = -{l + i-m^/A, for J =2; (5) 

fci = -fc2 + 77 +(1- (-1)^)0/4, for J =3. 

The degeneracies at the 1-pole and the 3-pole are related by symmetry, Eq. (||) . 

Let us now consider the special cases q = 1,2,3: The case q = l corresponds 
to (/i = and is trivial; the Bloch Hamiltonian has one eigenvalue, no crossing, 
and no charge transport. The case q = 2 (or equivalently, = tt) is already 
interesting. The Bloch Hamiltonian, Eq. (|^) reduces to the basic paradigm for 
Chern numbers — Berry spin 1 /2 Hamiltonian: 

{t2 COS ^2 ) 0-3 + 1 COS fci ) CTi + (^3 C0s(fci -I- ^2 - ?7) ) 0-2 , (6) 

with cTj the Pauli matrices. Since the matrix is traceless, levels cross when it 
vanishes. This gives Eq. (||) and is all of F(fcj), provided kj ^ 7r/2. At these 
special points F(fc2 = ^/2) is the two great circles ti = and ^3 = 0; similarly, 
F(fci = 7r/2) is the two great circles t2 — and = 0. If fci = fc2 = 77 = tt/2, 
then the whole unit sphere |t| = 1 belongs to F. 
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Now that the set of level crossings is known, we can describe the Chern 
numbers. By the general principles mentioned before, interesting Chern num- 
bers arise when the orbit in deformation space 7 traps level crossings. Let 7^- 
denote a small closed orbit around the j-th pole. For /c2 7^ 7r/2 these orbits trap 
level crossings and are such that the Chern number, Eq. (2), is well defined. 
The charge transport can be computed by a formula of 7]: 

Qili,k2,v) = Tsgn(cos/c2)sgn(sin(fc2 - 77)), 

Qh2,k2,v) = 0, (7) 

Q{l3,k2,v) = ±sgn(cosA:2)sgn(sin77). 

The overall sign depends on the orientation of 7^, and is opposite for the 
top/bottom bands. The Chern numbers change (discontinuously) on F so the 
direction of charge transport can be flipped by tuning k2 and ?/. 

For k2 = 7r/2, the Chern number (3(72, fe, rj) is not well defined since there 
are level crossings on the surface of integration. The total charge transport is a 
well defined integer if fc2 = it/ 2 is not an allowed channel, and is ill defined if 
this channel is allowed. For finite strips with periodic boundary conditions, odd 
strips have fc2 7^ 7i'/2 and the total charge transport is integral. Even strips, 
and also the infinite two dimensional lattice, include the k2 = 7r/2 channel and 
do not have a well defined (total) Chern number. For finite strips where the 
channel k2 — tt/2 is excluded, the total transport can be read off from Eq. (^. 
In particular, with maximal breaking of inversion symmetry, 77 — tt/2, an orbit 
of deformations 71 about the 1-pole, transports #{A;2 channels} charges in the 
ground state. The total charge transport is therefore a nonzero integer for any 
strip (where the number of fc2 channels is finite) and can be arbitrarily large. 
This shows that summation over the ^2 channels does not cancel in general: 
in this case, they add. In contrast, for an orbit of deformations 73 around the 
3-pole, the total charge transport is ±1 for all odd strips. This is because the 
allowed values of ^2 are equally spaced and then J2 sgn(cosfc2) = ±1. 

We see from this that: 1. The Harper model, Eq. (|l|), has nontrivial piezo- 
electric response. 2. For appropriate values of parameters and orbits, the charge 
transport is given by nonzero Chern integers. 3. The Chern numbers can sum to 
nonzero integers when summation over channels is taken. 4. Integer transport 
occurs also for arbitrarily small deformations 7. 

One may criticize the q = 2 example of piezoelectric response as being too 
special in that the deformations 7 that give charge transport are about points 
in parameter space where two hopping amplitudes vanish. This is a rare event, 
analogous to multicriticality. Can one have piezoelectric transport also if all 
hopping amplitudes remain positive? As we shall see, this happens for the 
Harper model with q = 3. The price we shall pay is that the analysis of the set 
of level crossing is more complicated and part of it relies on detailed numerical 
analysis. 

For q = 3, the model is described by a 3 x 3 matrix with the characteristic 
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polynomial 

- 3E ^ 2h{i, k, rf) = tl cos 3fci + (8) 
^2 cos3fc2 — ^3 COS (3(fci + fc2 — + itit2H cos(r/). 

Eq. (^) is a strong version of Chamber's relation: the coefficients of E are not 
only independent of fc, but also of t and rj. Therefore, the band edges are at 
extrema of h in the entire five-dimensional parameter space. The set of curves 
where the first gap closes for g = 3 is obtained when E = —1 and h = 1. 

The strategy we use to get hold of the degeneracy surface F is the following: 
At the 1-pole, Eq. (|) gives a two dimensional planar piece of F. (This two 
dimensional plane projects to a line in Fig. 1.) The line T(ti — 1,7] ~ 0) turns 
out to be a line of self-intersection of F. One two dimensional piece is given in 
Eq. (^). The intersecting two dimensional piece can be obtained as follows. Pick 
any point on F(ti ~ 1, 77 = 0) and expand h in powers of ^2 and require that 
/i = 1 to every order. If we now use k2 and t2 as the parametric representation 
of F we find: 

ki = sin(3fc2)i2/3 + sin(3A;2)i2 + ••■ ; 
tl = l-tl + cos(3k2)4 + ... ; (9) 
77 = -sin(3fc2)t2 - sin(6fc2)ti/2 -I- ... . 

In this parametric representation, ^2 is small and k2 is arbitrary. This gives us a 
thin strip of F, which intersects that of Eq. (||). We can extend this strip using 
the fact that the tangent plane to F is the kernel of the Hessian of h. In other 
words, with ^2 fixed, the curve of degeneracy, F(fc2), may be described as the 
solution of an ordinary differential equation: The velocity in parameter space is 
given by the kernel of the Hessian of h. Near the 3-pole, F(fc2) is given by 
and (|). 

Several curves describing degeneracies are shown in the figure. We chose ti 
and t2 as our coordinates on i restricted to the positive quadrant, < tj < 1. Let 
us denote by Fc the degeneracy surface with the poles excised. The tongue-like 
curves in Fig. 1 are Fc(fc2 — — 7r/18) and Fc(fc2 = — 7r/3). By the von Neumann- 
Wigner rule one expects these to be one dimensional curves, and indeed they are. 
For the orbit of deformation 7 shown in this figure, that is, a small circle centered 
at tl = 0.355,^2 = 0.446, the charge transport is (5(7, — tt/IS, 7r/9) = ±1 (the 
sign depends on the orientation on which 7 is traversed) . This gives an example 
where a Chern number is nonzero for a small orbit of deformations that lies 
entirely in the positive quadrant of hopping amplitudes. 

For g = 2, we have seen that Chern numbers for the infinite crystal included 
channels with ill defined Chern numbers. One may wonder if this is also the 
case for q = 3. The answer is no. For q — 3, all sufficiently small paths 7^ 
around the poles have well defined Chern numbers for all ^2 channels, and some 
of these are nonzero. It is easy to verify this for 77 = 7r/2 where one can check 
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that all of r is at the poles. More is true; T is in fact restricted to the poles 
for all 7r/6 < i] < tt/2. One way to see this, is by analysis at the vicinity of the 
point i = l/\/3(l, 1, 1), k = — 7r/18(l, 1), 77 = tt/6. One finds that rj attains its 
maximum value on Tc at this point. Further study shows that, in fact, for all 
77 7^ mod TT a sufficiently small orbit 7^ about the j-polc avoids T. The Chern 
numbers for these orbits are all well defined, and by numerical integration we 
found Qi'^j, k2,rj) = ±1, 0, =Fl for j = 1,2,3 pole respectively, for all k2 channels 
and T] ^ 0. 

In conclusion: we have described a method for analyzing the Chern numbers 
that arise in inversion asymmetric Harper models and have found explicit situ- 
ations with nonzero quantized piezoelectric response. In all these cases nonzero 
transport occurs for arbitrarily small orbits, and this can happen also when all 
the hopping amplitudes are positive. 
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Figure Caption: 

Curves on which the first gap closes for g = 3. The vertical axis is the 
asymmetry flux 77. The horizontal plane is the positive quadrant in the plane 
< ti,t2 < 1, tl+tl < 1. The hatched region is this positive quadrant at 
77 = 7r/9. The two vertical lines at ^2 = correspond to the gap closure in 
Eq. (|5|). The tongue-like curve is the line of gap closure at ^2 = — tt/IS. It links 
with a small circle 7 in the hatched plane, centered at {ti — .355, t2 — .446). 
A periodic deformation along 7 transports a unit of charge. The curve in the 
plane -q = Q is the line of gap closure for k2 = — 7r/3. 
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